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ABSTRACT 


Measurements of a 1230-km baseline were made during an eight-week period in the 
fall of 1987 using Trimble 4000SX single-frequency, five channel Global Positioning 
System (GPS) receivers. Twenty-eight days of carrier phase data were processed using 
correlated triple differences with fixed satellite orbits, the broadcast ephemerides, a 
modified Hopfield tropospheric model. and without ionospheric correction to determine 
the accuracies and precisions of the slope distance and baseline components. The data 
were processed in ever increasing observing sessions to determine the optimized obser- 
vation periods required to achieve various orders of geodetic accuracies. 

The accuracies of the slope distances were better than 1.0 ppm for any observing 
period. The day-to-day repeatabilities of the slope distance measurements were better 
than 1.0 ppm (2c) after 20 minutes of observations. Accuracies and repeatabilities (2c) 
of the baseline components were better than 10.0 ppm after 20 minutes of observations. 
The correlated triple difference results were on the order of previous GPS surveys that 
used higher resolution differencing or external timing aids. Discussions include the ef- 
fects of ephemeris, tropospheric and ionospheric errors, and dilution of precision. 

Observation periods and mean slope distance errors were reduced when observations 
started close to and included the infinite peak of the Position Dilution of Precision 
(PDOP). The smallest variances were associated with observations about the infinite 
PDOP peak. 
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I. INTRODUCTION 


A priori’ knowledge of the observation periods required to achieve specified orders 
of geodetic accuracies is important in planning efficient and productive geodetic surveys. 
While terrestrial survevs have specified field and processing procedures and standards to 
categorize the geodetic accuracies of surveys [Federal Geodetic Control Committee, 
1984], only recently have standards been proposed for surveys conducted with the Global 
Positioning Svstem (GPS) [Federal Geodetic Control Committee, 1986]. Among the 
proposed requirements are standards for the length of observing periods aud satellite 
geometry. 

Field studies bv Remondi [1984] and numerical simulations by Fell [1980], Langley 
et al. [1984], and Landau and Eissteller [1986] studied optimized observation periods, but 
for baselines less than 100 km. Cannon er al. [1985], Bock er al. [1984], Goad er al. 
[1985], Mader and Abell [1985], and Bertiger and Lichten [1987] conducted long baseline 
surveys, but did not study optimized observation periods. One of the objectives of this 
thesis 1s to fill the gap between the above studies, i.e., examine the optimized observation 
periods for a long baseline. 

The optimized times will be examined using the correlated triple difference carrier 
beat phase observable because of its insensitivitv to integer ambiguities and loss of lock 
of the GPS carrier by the receiver. Another objective of this thesis is to add to the body 
of triple difference accuracy testing following a recommendation by Remondi [1984, p. 
259]: “More testing 1s required to establish the full accuracy potential of the triple dif- 
rence method. ” 

GPS carrier phase and pseudorange measurements were made during an eight-week 
period in the fall of 1987 with Trumble 40005X single-frequency, 5 channel receivers. 
The long baseline is approximately 1250 km in length between the National Oceanic and 
Atmospheric Administration’s (NOAA) Western Regional Center located at Sand Point 
in Seattle, Washington, and the Naval Postgraduate School (NPS), Monterey, 
California. The baseline was determined by locating the positions of its ends by con- 
necting them by independent short baseline survevs from nearby Very Long Baseline 
Interferometry (VLBI) horizontal control points. The results of those surveys form the 


reference to which accuracy will be determined. 


Additionally, studies for repeatability were conducted following another recommen- 
dation by Remondi [1984, p. 263] to enhance the capabilities of GPS measurements. 
The recommendation was to perform extensive repeatability studies on non-varving 


baselines for verifying and improving the GPS modelling. 


ho 


ll. BASELINE DETERMINATION USING GPS 


A. INTRODUCTION 

The Department of Defense's (DOD) Global Positioning System is intended to 
provide accurate positioning and precise timing for navigation purposes by broadcasting 
codes superimposed on two radio carrier frequencies from satellites. The satellites are 
placed in a constellation so that at least four satellites are visible globally. The Precise 
Code (P code) will be limited to authorized DOD users. The Coarse Acquisition (C A) 
code provides real-time accuracies to about 100 m [Daker, 1986] and is available to 
anvone. 

Ihe codes provide their transmit times, satellite orbit and clock information, and 
information to enable anv receiver to lock onto other GPS satellites. The orbital infor- 
mation (ephemeris) provides the position of the satellite. The receiver measures the time 
delay between the receipt of the C/A code and its transmission time. The tirne delay can 
be transformed into an apparent slant range from the satellite's known position to de- 
termine the location of the receiver. Since it includes delays due to receiver clock errors 
and the effects of atmospheric refraction, the apparent slant range is referred to as the 
pseudorange. A minimum of four satellites are required to solve the system of range 
equations for the receiver's coordinates and clock errors. 

While the C/A code provides real-time location, it does not meet the accuracy re- 
quired for precise geodetic work. Nevertheless, GPS makes possible a higher resolution 
via the carrier signal. Though the carrier itself does not contain the orbital and timing 
information. Which would have to be supplied by some other means, it does offer a 


higher resolution because of its 19-cm wavelength. 


B. THE MONTEREY-SAND POINT BASELINE 
l. General 

The length of the Monterev-Sand Point baseline (Figure 1) was computed by 
differencing the World Geodetic System 1984 (WGS84) [Defense Mapping Agency, 
1987] Cartesian coordinates determined for Monterey and Sand Point by short baseline 
GPS surveys from known horizontal control points. The precision and agreement with 
terrestrial survey results of short baseline GPS measurements using single frequency, 
double difference solutions are well documented (e.g., [Remondi, 1984], [Goad and 
Remondi, 1984], [Bock er al., 1984)). 
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Figure 1. Monterey-Sand Point Baseline and environs: [nsets not drawn to scale. 


METEOROLOGY 
DEPARTMENT 


At Sand Point, on-site meteorological measurements were made near the middle 
of the observing session. For the Monterey antenna determination meteorological 
measurements were made every half hour and the mean of all the measurements was 
used in the processing. The 7rim640 solutions were obtained using uncorrelated double 


differences and estimating initial mteger ambiguities. The ambiguities were fixed to the 


integer values that produced the smallest residuals. A tropospheric factor was estimated 
along with the integer ambiguities and the baseline components in the least-squares 
processing. 

The horizontal control points used for the reference stations in determining the 
coordinates of the antennas were mobile Very Long Baseline Interferometry (VLBI) 
sites. Ihe NADS3 Cartesian coordinates for the VLBI sites were provided by the the 
Gravity, Astronomy and Space Geodesy Branch of the National Geodetic Survey (NGS) 
[Abell. 1987]. The NAD83 coordinates were determined in August 1987 from a global 
adjustment of VLBI surveys. Carter er al. [1985] described the NGS VLBI program. 
The Defense Mapping Agency Hydrographic’ Topographic Office (DMAHTC) validated 
the direct transformation of the VLBI Cartesian coordinates to WGS84 Cartesian coor- 
dinates [Kumar, 1988]. 

2. Mfonterey Coordinates 

The coordinates for the Monterey antenna location were determined bv aver- 
aging two surveys conducted on separate days from the VLBI site FT ORD NCMN 
1951. Table 1 lists the observing sessions used to determine the WGSS4 coordinates of 


the Monterey antenna. Table 2 lists the results of the GPS surveys. 


Table 1. MONTEREY ANTENNA LOCATION SURVEYS: From FT ORD 
٩٢17 ۱ ۰ 1981. PST = Pacific Standard Time. 


Number of 
Double BE EEN Slope 
Difference Rejected Distance (m) 


Observations 
09 16 87 1018 1210 996 7 12139.8125 02025 
09:18.87 0914 1120 1136 3 12139.8119 0.025 


The Monterey coordinates were computed from: 





Cury = Corn + 4, 


where: 
با‎ Monterey coordinate 
a Ft Ord coordinate 
A, Variance-weighted mean baseline component 


The uncertainties of the Monterey coordinates, oc,,,, , were computed from: 


Table 2. RESULTS OF MONTEREY ANTENNA LOCATION 
SURVEYS: WGSS4 Cartesian Coordinates (meters). 


09:16 S7 | -2697026.493 -10313.601 | 0.032] -2707340,093 0 A: 
09 18,87 | -2097026.493 -10313.602 | 0.037 


-2707340.095 A: 037 
Mean: -2707340.094 
09 16,87 | -4354393.309 | 0.010 | 917.693 | 0.048 ۱ 0.049 
09:18 87 | -4354393.309 | 0.010 | 917.689 | 0.049 


0 KE 
Mean: - 475.618 O. B 


09:16:87 3788077.778 0.009 -6337,391 | 0.0411 (۳۶ 0.042 
09 18/87 3788077.778 0.009 -6337.389 | 0.044] 3781740.389 0.045 


Mean: 3781740.388 | 0.044 | | 0.044 | 
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0.007 
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where oz, 1s the uncertainty of the Ft Ord coordinate. 

One month prior to the surveys originating from Ft Ord, two other survevs were 
conducted from the satellite Doppler horizontal control point NAVAL POST GRAD 
31965;DOPPLER. The Doppler station is approximately 300 m north of the Monterey 
antenna location (Figure 1). The double differenced GPS carrier phase solutions of the 
two Doppler-originating surveys yielded the following Monterey coordinates: 

۲ پ)ؤؤ‎ 777٦ه‎ 
Y  -4353475.654 m 
Z 3781740.264 m 


The baseline components had an uncertainty of +0.002 m, but the Doppler station has 
an uncertainty of +2 min each coordinate before transformation to WGS84. The results 
of the Doppler surveys were not used in determining the Monterey coordinates because 


of the large uncertainty in the Doppler station location. The three-dimensional positions 


of the Monterey antenna from the Doppler and the VLBI originating surveys agree to 
Deller tnan 0.5 m. 

The two days of pseudorange observations at the Monterey receiver were each 
subjected to a least-squares estimation and then averaged together to yield the Monterey 
coordinates: ۱ 

5009225095543 248 1.9 1 
Y | -4353466.807 40.6 m 
2 37/61/4110 70 


where the two-day standard deviations are given. The deviation of the pseudorange 
from the differenced carrier phase derived Monterey coordinates is expected because of 
the coarser resolution of the C/A code and because the pseudoranges are corrected for 
neither tropospheric nor ionospheric delays. 
3. Sand Point Coordinates 

The Sand Point antenna coordinates were determined by one 90-miñute GPS 
survey (Table 3) from the mobile VLBI site Aviation 2 which is 530 m distant 
(Figure 1). 


Table 3. RESULTS OF SAND POINT ANTENNA LOCATION 
SURVEY: From double difference carrier phase solutions. Distances 
and WGS84 Cartesian Coordinates are in meters. 


Aviation 2 Coordinate Sand Point Coordinate 


X= -2295347.760 0.017 | -408.361| 0.002 5777537562] 0.017 
Y= 9 0.028 330.200} 0.002 -3637699.228 0.028 
Z= 44 0.032 73,813} 0.003 4693482.777 0.032 

















The Sand Point antenna coordinates computed from the pseudorange data 


were: 
X  -2295749.623 m 
Y | -3637694.279 m 
2 1693488.472 m 


and standard deviations are not given because this is a single observing session and the 


program does not provide solution standard deviations. 


4. Mlonterey-Sand Point Baseline Components 
First the Monterev-Sand Point baseline components were computed bv sub- 
tracting the Sand Point Cartesian coordinate from the respective ۶ 
coordinate. The uncertainties in the baseline components were computed as the square 
root of the sum of the squares of the uncertainty of the Monterey and Sand Point co- 
ordinates. The slope distance was computed as the square root of the sum of the squares 


of the baseline components. The slope distance uncertainty was coniputed as: 
PANN لد‎ AA ) 
لد تس‎ ۳ OD Û 
where b is the slope distance. Using the information from Tables 2 and 3 gives 


Table 4, which is used as a standard to estimate accuracy. 


Table 4. STANDARD BASELINE DISTANCES AND 2-SIGMA VALUES. 


AX -411583.973 ۱ )0.19 ppm) 
Ay -715776.390 (0.16 ppm) 


AZ -911742.388 (0.12 ppm) 
Slope Distance 1230045.280 (0.09 ppm) 





C. THE ONE WAY CARRIER BEAT PHASE 

The development of the carrier beat phase technique and a model for its application 
are given by Remondi [1984]. The phase measurement is done by beating the received 
carrier with a local oscillator internal to the GPS receiver. The slant range from a GPS 
receiver to a GPS satellite can be modelled in terms of the time it takes the signal to 
travel or the number of cycles that occur between the satellite and the receiver. The 
range in cycles will consist of an integer and fractional number of cycles. When a GPS 
receiver locks onto the carrier signal, it can immediately measure the fractional part and 
begin counting subsequent integer cycles, but it cannot measure or account for the initial 
integer number of cycles that preceded the initial fractional part. These missing cycles 
bias subsequent measurements and are called the initial integer ambiguity biases. 

Ihe signal does not take a direct path to the receiver as it 1s refracted by the 
ionosphere and troposphere. Additional errors are caused by the satellite deviating from 
its predicted orbit, errors in the satellite clock, and errors in the receiver clock. Anti- 
cipating that the observables will be used in relative positioning, that a single frequency 


receiver will be used and ignoring other error sources, such as multipath, the one-way 


carrier beat phase, ¢,(r.s.i), observed at receiver, r, from satellite, s, at observation 


epoch, i, can be modelled to first-order as: 


Pr r,s,i)= PL) — HP Ar, 1) — D plr,s,i) 


$ tL, ded ma i p(r,s,i) J 
— t(r4) + A(r,s,1) 


(2.1) 


where: 
i Epoch identifier 
r Receiver identifier 
S Satellite identifier 
~ (i) Phase of received carrier signal 
(r,t) Phase of receiver generated carrier signal 
ye Transmitted frequency of carrier signal 
D(r,s,i) X Satellite-receiver slant range 
p(r,s.i) X lime rate of change of slant range 
C Speed of light 
ge Receiver generated carrier frequency 
zi) Tropospheric delay 


A(r,s,l) Initial integer ambiguity 

E . (i) 

«(0 7 $() = 577 

T "0+ 

2.9 = (9 + $ 

0 Mean clock offset for both receivers 


0) Clock drift between both receivers 


and parentheses do not indicate factors or functions, but simply enclose identifiers. 


Brackets indicate factors. 


D. DIFFERENCING THE ONE-WAY CARRIER BEAT PHASE 
l. Single Difference 
The single difference (SD) is formed by differencing carrier beat phase observa- 
bles from two receivers at the same observation epochs. Following Remondi [1984], 
taking the difference, expanding the ¢(r,/) terms of Equation (2.1), and expressing dif- 


ferences in the between-satellite and between-receiver phases as f6(/) gives: 


SDS) = PEST) o 
= n [(نیا)م - 2,5.80)م]‎ 
n 4 ] ۵)2.5,( — p(1,s.0 Jé(D 
- ]/ - AEM + fo) 
+ [902.50 + 150190 
E £D, — «0,01 4,51) — ACL, 


where the terms are the same as Equation (2.1). Single differencing reduces or eliminates 
satellite orbital and clock errors because they are common to both receivers. 
2. Double Difference 
A double difference (DD) is formed bv differencing single differences between a 


reference satellite, f, and another satellite at the same epoch: 
DD(h,s,) = SD(h,t) — SD(s,i) 


Because the differences at the same epoch are taken with the same reference 
satellite, the double differences for each epoch are correlated. The advantage of the 
double differences is that the clock dependent terms - £ó(i) and Lf, — A. JE(i) - are elim- 
inated. The significance of the removal of those terms is to reduce from nanoseconds 
to microseconds the timing accuracy required to achieve one cycle accuracy. The 
Trimble 4000SX achieves sub-microsecond accuracy by using the C/A code timing in- 
formation [Ashjaee, 1985]. 

3. Triple Difference 
A triple difference (7D) is formed by differencing the double differences for the 


same satellite pair at some integer number of succeeding epochs, /: 
ID(i) 2 DD(h,s,i - D) — DD(h, s,i) 


The advantage of the triple difference is that it eliminates all the time inde- 
pendent terms, namely the initial integer ambiguities, A(r,s,1), and becomes insensitive 
to the initial ambiguities and any cycle slips when the receiver loses lock. 

The disadvantages of the triple difference are: another level of correlation, loss 
of resolution, and reduced number of observations. The triple differences are already 


correlated with respect to satellite because of the underlving double differences, and are 
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further correlated with respect to time because consecutive triple differenced observa- 
tions will have the DD(h,s,i + 4 term in common. 

For short baselines, integer ambiguities can easily be resolved because unmod- 
elled errors are highly correlated between the two antenna sites and are mostly elimi- 
nated by the differencing. Algorithms can take advantage of the integer nature of the 
initial ambiguities and solve for them. At longer baselines the unmodelled errors are not 
as highly correlated and not eliminated by the differencing. These errors fold into the 
initial ambiguities, so that the ambiguities are no longer integers. In some cases the 
ambiguities cannot be resolved (e.g., Henson and Collier, 1986, Tables 1 and 2). 

Because of the advantages the triple difference offers for long baselines, [ use the 
triple difference scheme. While the triple difference can be decorrelated by forming a 
weight matrix [Remondi, 1984], only the correlated triple difference software was avail- 


able to me. 


E. ERROR EFFECTS 

For long baseline GPS surveys, the primary errors are satellite orbit errors, 
ionospheric and tropospheric delays [Remondi, 1984, and Beutler er al., 1986]. Orbit 
(ephemeris) errors are the result of the departure of the satellite from the broadcast 
ephemeris orbit. The ephemeris 1s a predicted orbit for the satellite. Orbit errors prop- 
agate directly into the baseline measurements when the GPS orbit coordinates are fixed 
in the processing [Hothem and Williams, 1985]. Orbit errors can be the donunant error 
source affecting the repeatability of long baseline measurements [Lichten and Border, 
1987). 

The magnitudes of baseline errors increase with increasing baseline length because 
of the increasing projection of the ephemeris error onto the baseline component [Fell, 
1980]. Estimates of the broadcast ephemeris error range from 25 m [Beutler er al., 1986] 
to 100 m [Wells er al., 1986]. The magnitude of the effect of the ephemeris error on 


baseline accuracy has been traditionally approximated [Beutler ert al., 1984, Equation 


(2. D]: 

E E 

vu (2.3) 
where: 


b X baseline length 


p Slant range (receiver to satellite) 


il 


E, error in baseline length 


E, error in slant range 


The slant range to a GPS satellite is about 20,000 km, which translates to a baseline er- 
ror ranging from 1 ppm to 4 ppm using Equation (2.3). 

It is expected tliat the ephemeris error for the Monterey-Sand Point baseline will be 
towards the lower end of the range because the ephemerides are uploaded prior to the 
satellites entering the Yuma Proving Ground [Russell and Schaibly, 1980] near the 
California border. The ephemeris linearization error specification is to 1 m per day 
[Wells ez al., 1986]. 

The ionosphere disperses the code (the group velocity) from the carrier phase (phase 
velocity) because the C/A code has a frequency of 1 MHz while the carrier signal upon 
which it is superimposed has a frequency of 1575 MHz. The effect 1s to increase the 
pseudoranges, but decrease the carrier phase derived ranges (Smith, 1987]. Field exper- 
iments [Beutler er a/., 1986] showed that 10nospheric dispersion shortens baselines on the 
order of a few tenths to perhaps 2 ppm. 

[onosplieric error 1s proportional to the Total Electron Content (TEC) aiong the 
signal path and the cosecant of the elevauon angle of the satellite [Smith, 1987]. Thus 
the error is greatest for low elevation angles and least at the zenith. Wells et ai. [1986] 
estimated that the range errors due to the ionosphere are from 150 m at the horizon to 
50 m at the zenith. 

[onospheric activity is a function of latitude, longitude, time of day, season, and 
sunspot activity. [onospheric activity increases towards the equator and towards the 
sunlit portions of the earth. Diurnally, it has a minimum near 0600 local time with a 
maximum around 1600 local. lonospheric activity increases with the peak in the sunspot 
cycle. The minimum in the current 11-vear sunspot cycle occurred in 1986. Upon these 
systematic characteristics sporadic ionospheric disturbances are superimposed.{Henson 
and Collier, 1986] 

The tropospheric error is proportional to the refractivity along the satellite-receiver 
path and proportional to the cosecant of the elevation angle (Martin, 1980]. The index 
of refraction, n, is the ratio of the speed of light in a vacuum to the speed of light in a 
particular medium, in this case the troposphere. Because the refractive index is a small 
fraction greater than 1.0, a more convenient unit to work with is refractivity, N, where 
N =(n—1)x 10°. The magnitude of the tropospheric biases range from 20 m for 10? 
elevation angles to 2.3 m at the zenith [Wells et al., 1986]. 
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Henson and Collier [1986] have shown that triple difference measurements are un- 
affected by path-dependent ionospheric bias errors, but by path-dependent gradient 
ionospheric errors, and Martin [1980] has estimated that the combined ionospheric and 


tropospheric gradient errors are on the order of meters per hour and are proportional 
to the cosecant of the elevation angle. 


III. DATA COLLECTION, PROCESSING, AND ANALYSIS 


A. EQUIPMENT CONFIGURATION 
1. Hardware 

A complete description of the Trimble 4000S X receiver is given by 6 
Navigation [1987a]. NPS operates three Trimble 4000SX GPS Surveyor receivers, of 
which two were used in this study. The 4000SX is capable of observing the C’A code, 
integrated Doppler and carrier beat phases of up to five satellites simultaneously. [ts 
ability to use the C/A code allows the receiver to be used as a stand alone navigation 
system which determines position using Doppler-smoothed pseudoranges and velocities 
[Ashjaee, 1985]. 

For precise relative positioning the 400USX can transmit its data through an 
RS232 port to a microcomputer for storage on floppy disk for post-processing. The 
4000SX's ability to use the C/A code allows it to decode the GPS navigation messages 
so that it can track satellites automatically and determine. Most importantly, it uses the 
C, A code in a time transfer mode to determine any offset and drift of its own clock and 
thus provide accurate time tags for the observations without an external atomic clock 
or synchronization with the receiver at the other end of the baseline. 

The receivers were left on continuously to allow unattended data collection. 

Multipath-resistant Trimble microstrip antennas were installed at both the 
Monterey and Sand Point locations. 

2. Software 

For relative positioning, the receiver is controlled from the microcomputer by 
Version D of Trimble's Datalogger program. The reference position (the geodetic coor- 
dinates of the antenna) and the particular options chosen must be entered into the re- 
ceiver via the receiver kevpad. The receivers were set to use the reference position height 
for point positioning when less than four satellites were available. 

Each observation session was initialized to log data when a minimum of four 
satellites were 15° above the antenna’s horizon. Five satellites were designated for each 
observing session. The software logs the observables and receiver clock parameters to 
a floppy disk every 15 seconds and the C/A code-determined antenna position and Po- 
sition Dilution of Precision (PDOP) every five minutes. The GPS navigation message 


is logged to a separate file at the beginning of the session. 
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B. SATELLITE OBSERVATION PLAN 
|l. Satellite Selection 
The same five satellite (SV) (6, 9, 11, 12, and 13) were used for the entire eight- 
week observation period. These five satellites were visible at both stations for over 100 


munutes, and of these five satellites four were visible for three hours (Figure 2). 


START DATE/TIME: 1987/10/17 15:21:15. DAY OF YEAR 290 TOH 573675. 
STOP DATE/TIME: 1987/10/17 18:16: 1. DAY OF YEAR 290 TOW 584160. 


DATA AVAILABLE: 


STATION: SAND POINT 
SV 
SV 
SV 


SV 


STATION: MONTEREY 
SV 
SV 
SV 
SV 
SV 





Figure 2. Satellite availability: Each dot represents 10 observations; each col- 


umn, 10 epochs. From Trim640 output. 


2. Position Dilution of Precision (PDOP) 

Trimble Navigation [1987b] recommends that observations include the time that 
the PDOP goes to infinity. The Federal Geodetic Control Committee [1986] notes that 
initial results from investigations indicate that best results may be achieved when the 
Geometric Dilution Of Precision (GDOP) is changing value during the observing ses- 
sion, and proposes that observing sessions start with a high GDOP and stop with a low 
GDOP. 

Position Dilution of Precision is a component of the GDOP. GDOP is a 
measure of how satellite geometry degrades point position accuracy (Jorgensen, 1984]. 
For computational ease in the navigation solution the GDOP is defined as the square 
root of the trace of the covariance matrix of the errors in position and time with the range 
errors set to one [Milliken and Zoller, 1980]. The role and definition of GDOP and 
PDOP in GPS point positioning were applied to GPS relative positioning, i.e., good 


PDOP would provide better accuracy than poor PDOP [King et al. 1985]. Good and 
poor PDOP are shown in Figure 3. 

Landau and Eissfeller (1986], using numerical simulations in which they assumed 
a full 18 satellite constellation and a receiver that could track those satellites that mini- 
mized GDOP, found that better accuracy for triple difference solution for a 68 km 
baseline generally corresponded to high GDOP. They used a more complete GDOP that 
included consideration of ionospheric, tropospheric and satellite position errors which 
are neglected in the conventional GDOP. 

The 4000SX does not record GDOP, but it does record PDOP every five min- 
utes. PDOP relates to GDOP as: GDOP? = PDOP? + TDOP?, where [DOP is the Time 
Dilution Of Precision, the error in the user clock bias multiplied by the velocity of light. 
The expected uncertainty in a GPS point positioning solution is a product of the PDOP 
and the expected slant range error. The difference between PDOP’s at Monterey and 
Sand Point remained less than 1.0 for the entire eight week observing period. The PDOP 
peaks at two times in an observing session (Figure 4) - 60 minutes and 150 minutes. 
The PDOP peaks occur near when the satellites lie in a common plane causing the sol- 
ution of the linearized range equations to diverge and the PDOP becomes infinite 


(Jorgensen, 198]. 


C. METEOROLOGICAL PARAMETERS 

Meteorological parameters are needed for the tropospheric correction model used in 
the processing software. On-site meteorological observations were not available for 
Sand Point. Instead observations from the Weather Service Office (WSO) of the 
Seattle- Tacoma International Airport (SEA-TAC) were used for Sand Point. SEA-TAC 
is approximately 27 km from Sand Point (Figure 1). The NPS Department of Meteor- 
ology routinely collects real-time hourly observations of sea-level pressure, temperature 
and relative humidity from the National Weather Service’s data network. Observations 
that pertained to GPS observing sessions were entered into a file that is accessed by the 
batch file building program. 

While the Federal Geodetic Control Committee [1986] has proposed the use of on- 
site meteorological parameters, researchers have had success using standard atmosphere 
parameters for satellite geodesy (Fell, 1976] or extrapolated meteorological data 
[Rathacher er al., 1986. 
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Figure 3. Poor and good PDOP: From King ert al. (1985, Fig. 3.2]. 
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Figure 4. PDOP versus time 


Hourly observations for Monterey were obtained from the NPS Department of 
Meteorology. The instruments were located approximately 300 m south of the 


Monterey antenna location (Figure 1). 


D. PROCESSING SOFTWARE 

A complete description of the Trimble supplied Zrimvec software can be found in 
Trimble Navigation [1987b]. The data was processed using the Trinvec Trim640 pro- 
gram, Revision AB. Trim640 limits processing to 700 epochs, so the first 700 epochs for 
each observing session are used. Sand Point was used as the reference station and its 
coordinates kept fixed in the least-squares processing. Sand Point was chosen as the 
reference station because four satellite availabilitv occurred later at Sand Point than at 
Monterey. This avoided having to load the not-in-common epochs from Monterey at 
each processing. 

Trim640 uses the C'A code derived positions obtained at the lowest PDOP for the 
initial estimates of the baseline components. Trimó40 culls the best C/A code position 
during the data loading. No ionospheric correction is provided by the software, and only 
the broadcast ephemeris can be used to compute fixed orbit satellite positions. [n the 
triple difference processing, the only parameters estimated by the least-squares process- 
ing are the baseline components, AX, AY and AZ. 

A modified Hopfield tropospheric model [Goad and Goodman, 1974] is used to 
correct the carrier phase delay caused by the troposphere. The correction is a function 
of the atmospheric refractivitv computed from surface meteorological values of pressure, 
temperature, and humidity, and the elevation angle of the satellite. Larger corrections 
are required for low elevation angles, as the signal travels a longer path through the 
troposphere. The model corrects for at least 90% of the tropospheric delay [Remondi, 
198]. Trim640 allows only one pressure, temperature and humidity entry for each site 


per session. 


E. PROCESSING PROCEDURES 
l. General 
To study optimized times, the data from each observing session Were segmented. 
Each successive segment contained 10 minutes more data than its predecessor. For ex- 
ample, for an entire observing session that started when four satellites were available and 
stopped when less than four were available provides 175 minutes of observations. The 
first segment will use the first 10 minutes of data , the sixth segment will use the first 60 


minutes, while the eighteenth will use all 175 minutes. For each segment, the entire 


processing was restarted from the data loading. Reloading the data for each segment 
takes considerably longer than using the Trim640 option to flag data for processing, but 
reloading was done so that the processing does not use a best C/A code position from 
later in the observing session. 

Convergence of the least-squares solution was achieved bv doing five iterations 
using every tenth triple difference formed from every tenth double difference, followed 
by five iterations decreasing the triple and double difference increments to five, and 
finally five iterations using all triple differences formed from all double differences. 

Trim640 rejects those observations whose residuals exceed a multiple of the 
mean residual. The multiple of the mean residual is known as the edit multiplier. 
Trim640 uses 3.5 as the default value for the edit multiplier. I used the default value for 
the initial processing. 

Any segment that had more than ten percent of its observations rejected or 
Whose solution slope distance standard deviation (o,) was greater than 10 m was re- 
processed. The reprocessing was identical to the initial processing except that before 
invoking the triple difference process the pseudoranges for both stations were subjected 
to separate least-squares adjustment. The pseudorange processing improves the CA 
code derived initial estimates for the baseline components and corrects the carrier beat 
phase time tags. The carrier beat phase time tags are computed from the C/A code 
times, and are earlier than the C, A code times. 

lf the pseudorange processing failed to lower the rejections to ten percent, the 
edit multiplier was increased until the rejections reached ten percent. A ten percent re- 
jection level was observed for a few sessions and always occurred within the first thirty 
minutes of observations. The data were transferred to the Naval Postgraduate School's 
IBM 3033 computer for analysis. The data were analysed and graphics produced using 
the APL-based GRAFSTAT program. 

To study the effects of reducing observation time, five case studies were under- 
taken in which the observation start times were changed for processing. Each case studv 
followed the processing procedures outlined above. 

2. Batch Processing 

Processing is performed in a batch mode. A batch file passes parameters to a 

template. Trimble supplies command files that tell the Trim640 to use the template pa- 


rameters in processing. 
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The batch file is built using the program Barbld (Appendix A). Barbld builds a 
batch file by providing the appropriate file names and start and stop 0 
computes the appropriate meteorological paraineters for each segment by locating the 
applicable weather observations from the weather observation file, interpolating values 
at the start time, computing running means from each hourly weather observation, then 
interpolating the running means to the stop time of each segment. Two mullibars were 
subtracted from the SEA-TAC sea-level pressure to compensate for the 20-m elevation 
above sea-level for the Sand Point antenna. 

Iniually, processing was done on an IBM XT with a math coprocessor anda 
hard disk. Processing the 18 segments of an observing session took ten hours of com- 
puting time. Later, processing was performed on an 80286 based microcomputer run- 
ning at 10 MHz with an 80287-8 math coprocessor that reduced the processing time to 
three hours. 

Two minor problems with Trim640 were discovered during the processing. 
First, large values in range differences were found when using the range differences 
rather than the pseudoranges to improve the C/A code positions. The data were for- 
warded to Trimble Navigation for evaluation and an error was found in their software. 
The error had no apparent effect upon carrier phase difference processing. Second, 
Trun640 1s incompatible with one or more of the AST Research, Inc. device drivers 
supplied with the 80286 microcomputer: ram disk, print spooler, and extended memory. 
Removing the drivers allowed Trim640 to execute normally. 

At the conclusion of the batch processing, the slope distance, the baseline 
components, their standard deviations (o,), the number of observations, the number of 
observations rejected, and the RMS cycle fits were extracted from the Trímó40 output 


file and collected into files that held the data for a particular segment for each case study. 


F. ANALYSIS PARAMETERS 
The statistical parameters that will be used to evaluate the results are the error, the 


sample mean and the standard deviation defined as [Davis er al ., 1981]: 
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where: 
۲ # Error for day J 
C, Measured component for day d 
۲ expected values 
E Mean error or bias 
c? Variance 
g Standard deviation of E, 


oz Standard deviation of E 


Accuracy describes the closeness between the measurements and the expected values 
[Davis er al., 1981]. The degree of accuracy is determined to the magnitude of the mean 
error (E). The repeatability of the measurements will be expressed in terms of 2e because 
it approximates the 95% confidence level for single-dimension measurements [Federal 
Geodetic Control Committee, 1986]. The slope distance and individual baseline com- 


ponents are one-dimensional measurements. 


G. DATA AVAILABILITY 

Observations were made simultaneously at Monterey and Sand Point for an eight- 
week period beginning 29 September 1987. Observations were made Tuesday through 
Saturday except the days after federal holidays. Forty observing sessions were con- 
ducted of which 28 were used in the analysis and are listed in Table 5. The remaining 
12 days of observations were not used in the analysis for various reasons, which are 
listed in Table 6. 

For brevity the observing days will be referred to by their Julian day. Times in 
Table 5 are given in Pacific Standard Time (PST) rather than Universal Coordinated 


Time (UTC) for ease in the later discussions on diurnal effects. 
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Table 5. DAYS USED IN THE DATA ANALYSIS 


s Number of 
Start Time End lime > a 
Date Julian Day Triple Difference 
MA... eo "5 Observations 


09 29 87 0833 
09 30 7 0830 
OUST ا‎ 0823 
10.0287 0822 
10:03 87 0817 
10 06 87 2 دراور)‎ 
100787 0801 
10.05. 57 0757 
10/09/87 0753 
10,10 87 283 0749 
10.14 $7 28 0732 
IOS 2 0721 
10/20/87 2 0707 
ST 0703 
1012257 0658 
10:24. 87 0651 
10:27:87 0638 
10 2857 0634 
10/29/87 0631 
702.87 XC 0613 
11.04:87 0605 
11 03 87 0601 

054] 

0536 

0528 

0524 

0455 

0438 





lable 6. OBSERVATION DAYS NOT USED IN THE ANALYSIS. 


Date Reason for not analvzing 


1015 87 No SEA-TAC weather observations 

10/16, 87 No Satellites at Monterev for first 10 minutes 
10/23/87 Slope distance c, ^ 10 m for first 10 minutes 
10/30/87 Disk error 

10551557 No SEA-TAC weather observations 

11/06/87 Unhealthv satellites 


11/07/87 Unhealthy satellites 
11/17/87 Unhealthy satellites 
11/18/87 Unhealthy satellites 
| ۱ Unhealthy satellites 
0۸ 7 Unhealthy satellites 
11/24/87 NO satellites at Monterey for first 20 minutes 





po 


Day 283 was processed with Monterey as the fixed reference station because the 
data set would not partition into 10-minute segments when Sand Point was used as the 
reference station. 

The days that were nussing one or two segments were excluded from the analysis, 
so that changes in the sample variances would not be due to unequal sample populations 
between the seginents. 

Reprocessing the first segment for 23 October failed to reduce the slope distance 
sigma to less than 10 m because that segment had only nine triple difference observa- 
tions because the receiver frequently lost lock on the satellites. That dav was not used 
in the analysis, even though its other segments had slope distance sigma’s less than 10 


m without reprocessing. 
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IV. RESULTS AND DISCUSSION 


A. GENERAL 


To study optimized times, five cases are studied: 


1 Process all data when four satellites at least 15° above horizon 


ty 


Begin processing when five satellites at least 15° above the horizon 


us 


Begin processing as ın 2, but delete fıfth satellite 


da 


Begin processing data 40 minutes later than in | 


Un 


Begin processing data 70 minutes later than in | 


Case l is essentially the processing of the full data set. Each of the other cases is a 
subset of Case 1. Trimó640 allows the user to designate at what time within the full ob- 
serving period that the data loading should begin. For Case 3, Trimó40's ability to flag 
data was used to exclude the fifth satellite (SV 12) from the processing. Each day of 
Table 5 was processed for each of the five cases. 

During the course of the discussions it will be necessary to distinguish between ob- 
servation periods and the time of the observations fixed with respect to satellite geom- 
etry. As the satellite geometrv (PDOP) begins with the availability of four satellites, the 
time of observations can be defined in terms of the Case | start time. Observation pe- 
riods are deternuned from the start time for each case. Times of observations will given 
as equivalent Case 1 times and is obtained by adding the case observation length to the 
case's time offset from Case 1 (Case 2 and 3’s offsets are 20 minutes; Case 4, JO minutes 


and Case 5, 70 minutes). 


B. ACCURACY 
1. Slope Distance 

The slope distance errors for Cases 1, 2, 4, and 5 are presented in Table 7. Then 
results show that accuracy to better than 1.0 ppm is achieved for any observation period, 
but that there are differences among the cases and with changing observation periods 
within each case. 

Cases 1, 2, and 4 exhibit similar behavior as the observation period increases - 
they become less positive (or more negative) as the observation period increases until 
they reach a minimum, then they reverse their trends and become less negative. Positive 


error indicates that the measured baseline is longer than the standard values, so that the 
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measured baselines are exhibiting an accordion effect as they shorten then lengthen as 
more observations are included in the solutions. While the minimum occurs after dif- 
ferent observing periods (140, 110, 90 minutes for Cases 1, 2, and 4 respectively), they 
occur at about the same absolute time with respect to the Case | start time (140, 130, 
and 130 minutes for Cases 1,2, and 4 respectively). As the start times occur later with 
each case, the errors for shorter observing periods become less positive and the ranges 
of the errors for each case decrease (the range of Case | errors is +36 to -18.4 cm; +7.8 
to -18.1 cm, for Case 2; and -10.3 to -15.1 cin, for Case 4). The error after the entire 
observation session decreases with later observing start times ( -14.8 cm, -9.1 cm, and 
-0. cm for Cases 1, 2, and d respectively). 

Case 5 behaves similarly to the previous cases except that its minimum occurs 
after only 10 minutes of observations and adding more observations causes the error to 
become positive. The error is largest after the entire observing session (38.4 cm). Cases 
|, 2, and 4 start their observations prior to the first PDOP peak that occurs at 60 min- 
utes (Figure 4) while Case 5 starts after the PDOP peak, so that starting the observa- 
tions close to the larger PDOP peak and including the PDOP peak observations can 
reduce tne error and the required observation period. The effect of the PDOP peak upon 
the mean slope distance error is readily apparent in comparing Cases 4 and 5. Case 4 
remains negative without the early positive error, and Case 5 remains positive as it lacks 
most of the observations from about the first PDOP peak. 

The first PDOP peak at 60 minutes differs from the second PDOP peak in that 
it has a higher value, is symmetric, and occurs farther from the four satellite observation 
start time than the second peak is from the four satellite observation stop time, ie., the 
second peak occurs with lower elevation angles than the first which implies larger 
tropospheric and ionospheric errors. 

The results for Case 3 (excluding the fifth satellite) slope distance errors are 
presented in Table 8. Case 3 was studied because the Case 2 results showed a decrease 
in the initial slope distance errors, and Case 3 was to study the effects of observing five 
satellites. The results show that using only four satellites makes a difference of only a 
few centimeters from the results using all available satellites. It should be noted that 
Case 3 uses only three satellites once SV 6 sets after 100 minutes of observations. The 


Case 3 minimum, -11.5 cm, was less negative than the Case 2 minimum of -18.1 cm. 
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Table 7. MEAN SLOPE DISTANCE ERROR FOR CASES 1, 2, 4, AND 5 


Observation Case 
Period 
[min] ppm cm ppm cm cm | ppm 


0:29 
ج0‎ 
0.18 
U12 
0.08 
0.04 
0.01 
-0.02 
-0.05 
-0.06 
-0.10 
-0.13 
-0.15 
-0.15 
-0.15 
04 
-0.13 
-0.12 
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Table 8. CASE 3 MEAN ERRORS 


Observation 
Period 
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2. Baseline Components 

The AX and AY accuracies are better than 10.0 ppm for all observing periods 
while AZ accuracy is better than 1.0 ppm for all observing periods. The accuracv of the 
baseline components is expected to be less than the accuracy of the baseline because the 
baseline errors are mostly perpendicular to the baseline itself [Remondi, 1954]. The 
baseline component results (Tables 9, 10, and 11) show that the AX and AY errors are 
ereater than the baseline components while the AZ errors are about the same order of 
magnitude as the slope distance errors. The AX and AY errors are negatively correlated 
which 1s the result of the correlations of the triple differences and the semi-circular tracks 
of the satellites (Figure 5). Case 3 (Table S) shows little difference from Case 2 in the 
AX error, and a more negative AY error is offset by a less negative AZ error. 

The smallest mean errors for the baseline components are found at various ob- 
serving periods. Zero mean error for all the baseline components is achieved with fewer 
observations as the observing start times occur closer to and before the larger PDOP 
peak. 


Table 9. MEAN AX ERROR FOR CASES 1, 2. 4, AND 5 


Observation 
Period 
S 173 
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Table 10. MEAN AY ERROR FOR CASES 1, 2, 4, AND 5 


Observation Case l Case 2 Case d E 
Period 
[min] ss ppm ppm 
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Table 11. MEAN AZ ERROR FOR CASES 1, 2, 4, AND 5 


Observation 
Period 
= [= [we | [om | جر‎ [me 
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Figure 5. Sky Plots of satellite tracks for Monterey: Elevation angles are dotted 


concentric circles. Zenith is at the center. 


C. REPEATABILITY 
I. Slope Distance 

The day-to-day repeatabilities, represented as the 2o level, for Cases 1, 2, 4, and 
5 are presented in Table 12, and for Case 3, in Table 13. All the cases achieve 1.0 ppm 
repeatability for anv observing period except Case | which requires 20 minutes of ob- 
servations. Repeatability eventually reaches better than 0.5 ppm after 60 nunutes of 
observations for any case. 

The minimum 2o levels for all cases are reached at the 80 to 90 minute time of 
observation, which is about 30 minutes after the larger PDOP peak. A slight increase 
in the 2a level is centered about the 120 to 130 minute time of observation for all cases 
which is near the PDOP minimum, after which the 2c level decreases slightly as obser- 
vations from the second PDOP peak are included in the solutions. Case J had the 


narrowest range of 20 values, and Case 1 had the widest range of 2c values. 


Table 12 SLOPE DISTANCE 2-SIGMA VALUES FOR CASES 1, 2, — AND 5 
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Period 
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Table 13. CASE 3: 2-SIGMA VALUES 


Slope ; 


e | سم‎ a چتھتھ نچ‎ ER 










Observation 
Period 
[mun] 


















62.4 76.6 
59. ۱ 80.8 
57-8 76.6 
58.0 د څ‎ 
58.8 9 4 
54.0 E TS 97.0 
53.0 279 95.8 
53.6 2.79 $9.0 
56.4 2.1 83.4 
58.0 2:32 72 
59.2 1.94 TES 
59.0 98 08.2 
59.0 1.71 65.0 
58.8 iS 65.4 
58.8 1.79 64.8 
5چ‎ 1.91 63.8 









2. Baseline Components 


All the baseline components have repeatabilities better than 10.0 ppm for any 
observing period except for the Case | AX, which required 20 minutes of observations 
(Tables 13, 14, 15, and 16). 

[t 1s interesting to note that while the AX 2c values for the first segment of 
Cases 2, 3, 4, and 5 are less than Case l's first segment, the final Case 1 2o value is less 
than the final segment of any other Case. The final 2e values for Cases 2, 3, 4, and 5 
are greater than the Case | 2o values after an equivalent number of observations. 

The minimum 2c levels for AX and AY occur at or near the end of the observing 
sessions for Cases 1,2, 4 and 5. For Case 3 the minimum 26 levels occur after 130 min- 
utes of observations. The minimum AZ 2c level occurs after various observation peri- 
ods, but generallv in the vicinity of 130 to 150 minute observation time, which 1s between 
the PDOP minimum and the second PDOP peak. Case 3 behaves in an opposite fashion 
from the other cases in that its minimum AX and AY 2o levels occur at the 150-minute 


time of observation while its AZ minimum occurs at the end of the observation period. 
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Table 14. AX 2-SIGMA AA FOR CASES 1, 2, 4, AND 5 


Observation 
Period 
[mun] 





Table 15. AY 2-SIGMA VALUES FOR CASES 1, 2,4, AND 5 

Observation 
Period 
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Table 16. AZ 2-SIGMA VALUES FOR CASES 1, 2,4, AND § 


Observation [Case 
me 
[min] ppm cm ppm 
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3. Standard Deviation of the Mean 

The repeatabilitv values can be used to estimate the standard deviations of the 
mean errors given in Tables 7 through 11 by using Equation (5.2). The values of Tables 
12 through 16 should be divided by 28 (where 28 is the sample population) to compute 
the standard deviations of the means (at the 2c level). Generally, the repeatabilities were 
about five times the magnitudes of the mean errors; therefore, the uncertainties of the 
mean errors are on the order of the mean errors themselves. Allowing for the 0.1 ppm 
uncertainty in the baseline and the baseline components and for the possible standard 
deviation of the means, accuracies to better than 1.0 ppm for the slope distances and 


10.0 ppm for the baseline components remain valid. 


D. ERROR EFFECTS 
1. 7-Day Means 
Because of the observations were made over a long period of time, Case 1 was 
subdivided into four groups comprised of seven consecutive observation days to study 


trends in the slope distance error to identify the contribution of various error sources to 
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Optimizing observing times. The results for the slope distance error are presented in 
Table 17. 


Table 17. CASE I ERROR: 7-DAY MEANS 


Observation 
Period 
[min] ppm 
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3.4 
MS 





The seven-day mean slope distance errors remain below 1.0 ppm for all groups 
for any observation period. The differences between the groups are mainlv in the 
predominance of negative errors in Groups 2 and 4 while Groups l and 4 have pre- 
dominantly positive errors. 

To account for the change in the characters of the groups, the change in the 
major sources of error are examined. 

2. Ephemeris Errors 

Examining the AODE (Age of Data Ephemeris) for all satellites during the 
eight-week observing period, showed that the oldest AODE for a single satellite was nine 
hours, but most were less than five hours. Averaging the AODE for all satellites for 
each observing day showed a range of seven hours during the eight weeks of observa- 
tions (Figure 6). Assuming 12 hours is the largest AODE, the maximum ephemeris 
linearization error would be 50 cm. Using Equation (2.3) results in 0.025 ppm baseline 


relative error. 
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Figure 6. Mean age of data (ephemeris) 


The changes in the seven-day mean AODE’s do not correspond well to the 
changes in the groupings, especially as the largest change in AODE from group 3 to 4 
does not correspond to a similar change in the group 3 to 4 mean slope distance errors. 
This is not surprising because while the change is relatively large, the magnitude of the 
orbital errors 1s expected to be small. 

The ephemeris error will appear as a bias during any one observing session. Dy 
averaging over eight weeks some of the errors will cancel and appear as variance. Be- 
cause of the short observing sessions, little variation in the ephemeris error 1s expected 
during any single sessions 

3. 10105016116 5 

All the observing sessions completed well before 1600 local standard time and 
approached the normally ionospheric activity minimum of 0600 (Figure 7). The near 
north-south orientation of the Monterev-Sand Point baseline places both stations in the 
same time zone, so that ionospheric errors between the two stations will be correlated 


to some degree and reduced in the differencing. 
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Figure 7. Difference between observation end time and 0600 PST 
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Figure 8. Electron fluence: — Units electrons — cm? —day-! —sr-!. 
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The weekly Preliminary Report and Forecast of Solar Geophysical Data |Space 
Environment Services Center, 1987] for the eight weeks of the observations described the 
solar and geomagnetic activity as generally quiet or low. 

Though the observing period occurred during a general lull in ionospheric ac- 
tivity, there are two ionospheric effects that can be easily examined: the diurnal effect 
as the observing periods occur closer to the diurnal ionospheric activitv minimum and 
sporadic effects from electron fluxes. The change in ionospheric range error for the 
eight-week observing period was computed from Henson and Collier [1986, Equation 
(3)} using the difference in Total Electron Content (TEC) determined from Spilker [1980, 
Figure 1-11] and Henson and Collier [1986, Figure 2] for the inean observation times for 
the first and last groups. The maximum change in range error, because of diurnal 
ionospheric changes from the first to the last week, is approximateiv 1.5 m which by 
Equation (2.3) is equivalent to a change in the relative baseline error of only 0.075 ppm. 

While the 0.075 ppm trend is not discernible because of the 0.1 ppm uncertainty 
in the baseline and the larger uncertainties in the measurements, the mean number of 
observations (Table 5) increased with each later group. The mean number of triple 
difference observations available were : 1519 #45, 1522 +22, 1601 +124, and 1726 +174 
for the first through the last groups respectively. The increased number of observations 
could be attributed to stronger signal-to-noise ratio with the reduced ionospheric 
dispersion as the observations approach the diurnal ionospheric activitv mununum. 

The group to group changes did not show an apparent trend, but alternated 
between predominantly negative and positive mean slope distance errors indicating a 
sporadic effect, or complex interaction between the error sources. To examine sporadic 
ionospheric effects, daily Electron Fluence (EF) from the weekly Preliminary Report and 
Forecast of Solar Geophysical Data [Space Environmental Services Center, 1987], and 
seven-day means of the EF were computed (Figure 8). Electron fluence 
( electrons — cm? —day^ —sr ) is the daily average of electron flux with energies 
greater than 2 Mev as measured by the GOES-7 satellite. High values of the seven-day 
mean EF corresponded to predominantly positive seven-day mean slope distance errors 
(Groups | and 3), and low mean EF corresponded to negative mean slope distance error. 

4. Tropospheric Errors 

Refractivitv, N, was computed using the meteorological parameters from the 

time closest to the middle of each observation session and Remondi (1984, equations 


2.30 and 2.31]. Possible indicators of the influence of the troposphere upon the baseline 


3] 


solution are the difference in refractivity, AV, between the baseline stations (Figure 9) 
and the mean N of both stations (Figure 10). The difference can indicate the level of 
correlation of the tropospheric errors. Correlated tropospheric errors are reduced in 
differencing. The mean N (N) will indicate the magnitude of the refraction of the carrier 
signal to both sites. 

The changes in the difference in the seven-day mean NV (Figure 9) did not have 
an apparent correspondence to the changes in the seven-day group mean slope distance 
errors. Instead N shows a pattern similar to the group slope distance errors 
(Figure 10). High mean refractivity corresponded to positive mean slope distance errors 
(Groups I and 3), and low mean refractivity corresponded with negative mean slope 


distance errors. The difference between the high and the low mean N are small and do 


not readily account for the magnitude of the changes. 





Figure 9. Difference in refractivity between Monterey and Sand Point 
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Figure 10. Mean Monterey-Sand Point refractivity 
5. 7-Day Repeatability 
The Group repeatabilities could not be given a general characterization as the 
Group means were because the was no consistent pattern to the Group-to-Group 
changes in their repeatabilities. For completeness, the 2¢ values for the various error 


sources are presented in Table 18, and the Group 2c values are presented in Table 19. 


Table 18. GROUP 2-SIGMA VALUES FOR THE ERROR SOURCES 
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Table 19. CASE 1; 2-SIGMA VALUES, 7-DAY GROUPINGS 


Observation 
Period 
min] em | ppm ppm | cm 


1.89 
1.01 
0.73 
0.53 
0.43 
0.37 
0.36 
0.39 
0.39 
0.41 
0.43 
0.45 
0.47 
0.49 
0.49 
0.48 
0.48 
0.47 
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E. EFFECTS OF THE C/A CODE 

The best C/A code position for Monterey was extracted from each observing session 
for comparison with the carrier phase results. The coordinates of Sand Point were sub- 
tracted from the 28-day mean of the Monterey C/A code coordinates to form the 
baseline components. The slope distances were computed from the square root of the 
sum of the squares of the baseline components. The mean errors of the C/A code de- 
rived baseline Were on the order of the carrier phase mean errors, but the 2c levels are 
several times greater than the carrier phase results (Table 20). These results are for one 
instantaneous position determination and would improve with time averaging within the 
observing period or Doppler-smoothing. 

The C/A code has its most significant impact upon the differenced carrier phase 
solutions when there are few carrier phase observations. Table 21 shows that reproc- 
essing the carrier phase observations was required for several segments because either 
more than ten percent of the observations were rejected or the slope distance solution 
sigma was greater than 10 m, and that all the reprocessed segments were confined to the 


first two segments. 
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Table 20. BEST C/A CODE RESULTS AND ERRORS 


























Component Distance (m) Relative Error (ppm) 
1 -411587,334 "I. RU 
OT -715776.997 SS TAS 
AZ -911759,358 SS TT 
Slope Distance 1230044.511 DO +t 72 





Table 21. CASE 1: DATA SEGMENTS BEFORE AND AFTER REPROCESS- 
ING 


| | Befoe | After 
Total 
Segment Slope ۱ ۱ Slope Ed 
Obs SEHE Distance | ی‎ Distance E. 
Obs Multiplie 
(m) (m) 


12500592192 
1230268.988 
1230047.276 
1230227.492 
۱ (9 8 
1230048.636 
1230044.230 
1230042.920 






































1230045.846 
1230045.294 
1230047.020 
1230045. 717 
1230045.915 
1230045.326 
1 6 
1230045.972 
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281 2 1230034.198 0 1230045.082 کیہ‎ 
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Case 5 
287 P. 43 1230043.418 
308 24 1230043.047 
314 | 1230051.620 


The carrier phase solutions Were improved by the least-squares processing of the 
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pseudoranges to obtain better initial antenna position estimates and also to correct the 
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carrier phase time-tags. All the Case | reprocessed segments were improved bv forcing 
the processing program to accept more carrier phase observations by increasing the edit 
multiplier. 

Because of the success of improving the solutions by the pseudorange pre- 
processing, all segments for three davs (276, 287, and 294) were reprocessed to see if their 
already acceptable solutions would improve. The solutions either did not change or 


became slightly less accurate. 


F. COMPARISON WITH PREVIOUS STUDIES 

Remondi [1984] concluded that 100 ppm accuracy was achievable in 30 minutes and 
10 ppm, in one hour - regardless of baseline length - based upon single-frequency, single 
difference solutions using precise ephemerides over baselines less than 100 km. To study 
the required observation times for various accuracies, he partitioned his data into 
15-minute observation spans. While decorrelated triple differences would be capable of 
performing relative geodesy at 1 ppm level, correlated triple differences may provide 5-10 
ppm or better [Remondi, 1985]. 

Some long baseline survevs include: 


Bock er al. [1984] achieved 2 ppm accuracy in closure of a transcontinental net aided 
by external atomic clocks using 10 hours of single frequency, single difference obser- 
vations. 


Cannon et al. [1985] achieved 0.1 to 0.7 ppm repeatability using single frequency 
non-differenced network solutions without ionospheric correction for 1700-km 
baselines between two California stations and Calgarv, Alberta, Canada over a four- 
dav period. The range in repeatability was attributed to satellite clock errors. They 
also found that the largest errors occurred in the AX component. 


Mader and Abell [1985] found that single frequency long baseline GPS results using 
single differenced carrier phase observables agreed to one ppm with VLBI measured 
baseline lengths. Their 2-day repeatability was 0.24 ppm. 


Goad er al [1985] measured a 1302-km baseline between California and Texas to better 
than 1.0 ppm: compared to the VLBI measured distance, and their two-day repeat- 
ability was better than 0.5 ppm. 


Bertiger and Lichten [1987] achieved 6 parts per billion repeatability over a 1314-km 
east-west baseline over separate four and seven-day periods using dual-frequency re- 
ceivers with a fiducial network for orbit determination, and some water vapor 
radiometers. 


Except for the Bertiger and Lichten [1987] results, the Case 1 results are competitive 
with the above studies despite the inherent low resolution and the added burden of the 
correlations of the correlated triple differences. A possibility for the optimum results of 


Case | is that a good solution is insured by setting both receivers to their known WGS84 
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coordinates. For real-time surveys, one end of the baseline will usually not be known 
to one meter accuracy. 

During the eight-week observation period. the third 4000SX was delivered to NPS. 
The antenna was installed on the same rooftop, but approximately 35 m to the nerth 
of the Monterey antenna. Two days of data using Monterey as- the reference station 
were used to fix the location of the new antenna. On Day 294, the new receiver's coor- 
dinates were offset 37 m, equivalent to 30 ppm baseline error. The difference between 
the offset antenna and Monterey solutions for the baseline to Sand Point is 0.7 ppm for 


the first seginent (Figure 11). 


0.6 


0.4 


PARTS PER MILLION 


0.2 


40 80 120 160 
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Figure 11. Results of station 2 offset 


For the remaining segments the difference is about 0.05 ppm. The first segment will 
be more affected by receiver offset because the small number of carrier phase observa- 
tions will be sensitive to the initial position estimate provided by the C/A code. The C/A 
code, in turn, will be sensitive to the receiver coordinates because the receiver will com- 


pute the slant range to the satellites using those coordinates. 


G. MISSING EPOCHS 
Missing epochs prevented a sixth case study in which the number of observations 


are reduced by flagging the data. The sixth case study was to have studied the ellects 
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FRACTION OF AVAILABLE OBSERVATIONS 
9.58 0.60 


0.52 


0.48 
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Figure 12. Fraction of available triple difference observations 


of reducing the number of observations without changing the geometry. Several days 
of data were processed in which only every other epoch was to be used in the processing. 
An hour of data was needed to obtain enough triple difference observations to consist- 
ently have less than 10% of the triple difference observations rejected and slope distance 
o's less than 10 m. This inability to form a reliable solution indicated a lack of carrier 
phase observations, so the number of observations were counted. 

The mean number of observations (accumulative) available for each segment, as a 
fraction of the maximum number of triple difference observations remained steady until 
120 minutes when it increased (Figure 12). The maximum number of triple difference 
observations is: Mr — (n; —l)(n, —1) where: 

JIr, maximum number of triple difference observations 


n number of observations epochs 


I 


n number of observed satellites 


$ 


Solutions using only 10 minutes of data at a time were examined to investigate the 
rise in available observations after 120 minutes. A curious phenomenon was discovered 


that while SV 6 was being tracked at Sand Point between five and ten of every 40 epochs 
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Were lost (Figure 13). As soon as SV 6 was no longer being tracked at Sand Point. no 


complete loss of epochs was observed (Figure 14). 


START DATE/TIME: 1987/10/21 16:43:07. DAY OF YEAR 294 TOW 319380. 
STOP OATE/TIME: 1987/10/21 16:52:60. DAY OF YEAR 294 TCH 319980. 
DATA AVAILABLE : 
STATION: SAND POINT 

SV 6 | 

SV 9 | 

5۷11 | 

5۷12 | 

دا | 5۷15 

STATION: MONTEREY 

SV 6 | 

SV 9 | 

5۷11 ا‎ TIME MEER ME S EE EE ie ehe 
E TL Lucus oer بی و ا‎ a و‎ 
5۷135 | .. 


Figure 13. Missing epochs while tracking SV 6: Each dot represents one obser- 
vation: each column, one epoch. Sand Point is missing epochs 4. 8, 14, 
18, 26, 30, 34, and 40. From Trini640 output. 


START OATE/TIME: 1987/10/21 17:03:04, DAY OF YEAR 320580. 
STOP DATE/TIME: 1987/10/21 17:12:60. OAY OF YEAR 321180. 
OATA AVAILABLE 

STATION: SAND POINT 

SY 9 | 

SV11 | 

ا 5۷12 

٦‏ د5713 


ا 6 5۷ 

STATION: MONTEREY 
09 de 

5۷11 | 

5۷12 | 

SVS aE cars 

SV 6 | 





Figure 14. Continuous tracking without SV 6.: No missing epochs at Sand Point. 


The phenomenon is not multipath interference because it affects all satellite at each 
epoch. A simultaneous outage would suggest interference from some other emitter at 
or near 1540 MHz, if the outages continued throughout the observation session rather 


than being limited to when SV 6 was being tracked. No physical cause for the outages 
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is readilv available, and the data have been forwarded to Trimble Navigation Ltd for 


further evaluation. 


H. DILUTION OF PRECISION AND RANGE ERRORS 

The case studies showed that the accuracy and precision of the slope distance were 
improved by using observations associated with the infinite PDOP peak rather than ei- 
ther the PDOP minimum or the secondary PDOP rise near the end of the observing 
sessions, and that improvements in the slope distance accuracy could be achieved using 
fewer observations when observations are made through the infinite PDOP peak. Pre- 
cision and accuracy of the solutions are dependent on geometry (fixed with respect to 
Case | time) and error sources rather than the number of observations (fixed with re- 
spect to length of observing period). 

The mean slope distance errors were governed by the ionospheric and tropospheric 
errors opposing each other. Early and late observations had positive mean errors while 
the errors about the infinite PDOP peak were negative. Tropospheric errors are partially 
corrected by the modified Hopfield model while the ionospheric errors remain uncor- 
rected, and as the Satellites reach the infinite PDOP peak the full impact of the 
lonospheric errors project upon the baseline. 

The ròle reversal of the PDOP peak from point positioning to relative positioning 
requires some investigation. Because PDOP reflects the strength of the satellite geom- 
etry for a single station, the investigation should begin with the satellite-baseline (two 
station) geometry by examining a simple projection of errors upon a baseline using a 
single satellite whose track is directly over the baseline (Figure 15). 


From the law of cosines: 
2 2 2E 
br = وم + ۱م‎ - 2۵۱۵2 7 (6.1) 


where: 
b, true baseline 
p, range from station | to the satellite 
p, range from station 2 to the satellite 


y angle between p, and p, 
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Figure 15. Relative positioning geometry: Not drawn to scale. 


Similarly we define a measured baseline, by: 
b3, — (pı + Eu)’ + (m + Epa) — جډیه)ل رت ې م)2‎ Ep) cos y (6.2) 


Where: 
E, error in p, 


۰۰۱۲ 


Subtracting the square root of (6.1) from the square root of (6.2), and letting 


cosy = |, gives: 
bez AE? ApAE 
لاک‎ To eee aig p. = (6.3) 
br سا‎ 


where AE, 2 Ej — E,, and Ap = p,— P, - lrom (6.3), Ap acts as an amplification factor 
for the difference in the range errors ( AE, ). As a satellite approaches the mid-point 


of a baseline, Ap 2 0 and b,» » AE, 
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Similar arguments can be extended to point positioning using two satellites at one 
epoch or two epochs of one satellite by inverting the geometry of Figure I[5, but the 
reduction in range errors would apply to point positioning as Well. 


The effects of y have been ignored because y is small. From the law of sines : 


b sin 0 


E (6.4) 


(7 = 
where O =e+v. e is the elevation angle to the satellite and v is the vertical angle be- 
tween the slope distance and the local horizon. Expanding sin? as 
sin e cos v + cos e sin y , and ( sin e cos v) > > ( cos e sin vy) because cos vx] ( v 1s -5? for 
the 1230 km baseline ), then sin 0 can be approximated by sine. Approximating b by 
1230 km and p bv 20000 km. then, by the law of sines, y is approximately 1? when e is 
15°; bv the law of cosines, y is approximately 2.5° when the satellite is at the baseline 
nudpoint. 

With cos 1°+0.9998 and cos 2.5°%0.9991, a very insignificant change, most of the 
deviation of the measured baseline from the true baseline must come from the error dif- 
ference factor, AE,. Most importantly having satellites at high elevation angles with re- 
spect to the ternuni of the baseline is not detrimental to relative positioning as high 
elevation angles are for point positioning. 

A second area of investigation 1s the path dependent errors. The ionospheric and 
tropospheric errors, approximated by cosecant functions of the elevation angle, are 
minimized at high elevauon angles. Simultaneous high elevation angles for both stations 
is another by-product of high PDOP. 

Finally, the effects of differencing must be explored. Qualitatively, when the satel- 
lites are bunched close together at high PDOP (Figure 3), some of the correlation in 
errors that was lost in lengthening the baseline is restored. The change in tropospheric 


and ionospheric errors with time will also be minimized because: 


d csc e — cose de 
di ww dt i 
and high elevation angles minimize Equation (6.5). Such time derivative minimization 
of errors, will greatly affect the triple difference because the triple difference is itself a 
time difference. 
The sum of the elevation angles for all satellites at both stations (Table 22) reaches 


a maximum after 70 minutes of observations, which is ten minutes later than the PDOP 
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peak. The minimum variance is achieved in the 80 to 90 minute Case 1 time. The run- 
ning mean of the sum of the elevation angle peaks at the 100 minute observing time 
While the running mean elevation angle reaches a local maximum in the 80 to 100 minute 
range. The running means provide a better estimator of the behavior of the accuracy 


and precision because of the accumulation of observations. 


Table 22. ELEVATION ANGLES 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Following a recommendation by Remondi [1984] that more testing of the accuracies 
of triple difference carrier phase measurements was needed, I studied the optimized GPS 
observation times required to achieve geodetic accuracies by partitioning observation 
periods for a 1230-km baseline into ten-munute segments and changing the length and 
starting time of observations. Accuracy was determined by comparing measured 
baseline values with reference values obtained by locating the ends of the baselines from 
very precise VLBI horizontal control points using double difference GPS carrier phase 
measurements over short baselines. 

The 1230-km baseline was directly measured for 28 days over a period of eight weeks 
using Trimble JOOUSX receivers. The Trimble supplied Trimvec software was used to 
process the carrier phase measurements. Using broadcast ephemerides in a fixed orbit, 
triple difference carrier phase solution with no ionospheric corrections and a 
tropospheric delay model that used only surface meteorological data, 1.0 part per million 
(ppm) accuracy in the slope distance was achieved for any observing period with a day- 
to-dav repeatability better than 1.0 ppm (2e). The AX, AY, and AZ (the components 
of tlie baseline parallel to the WGS84 axes) achieved accuracies better than 10.0 ppm for 
any observing period. AY and AZ repeatabilities were better than 10.0 ppm for any 
observing period, while AX required 30 minutes of observations to reach 10.0 ppm. 

[ had not expected to achieve 1.0 ppm accuracy in the slope distance measurements 
because Remondi [1984] had suggested that 1.0 ppm required dual frequency measure- 
ments using highlv accurate orbit information and water vapor measurements, and my 
solutions are further burdened by the correlations of the triple differences. The 1.0 ppm 
accuracies were on the order of the uncorrelated single frequency results of previous long 
baseline GPS surveys: Cannon er al. [1985], Bock er al. [1984], Goad er al. [1985], and 
Mader and Abell [1985]. 

The unexpected 1.0 ppm accuracy can be attributed to low ionospheric activity be- 
cause of the orientation of the baseline, the time of year, and the minimum in the 
sunspot cycle and solar activity during the observing period. Ephemerides errors are 
considered low because of the low age of the ephemerides. Those ameliorating factors 


must be considered before applying the results of these case studies to planning surveys 
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on different baselines at different times of day and year. lonospheric errors were on the 
order of the tropospheric errors. The superposition of the opposing ionospheric (short- 
ening) and tropospheric (lengthening) errors reduced mean slope distaiice errors. 

Slope distance errors and observation periods were reduced when GPS observations 
started near the infinite, svminetric PDOP peak. The reduction of errors about the 
PDOP peak can be attributed to the siniultaneous minimization oÉ range errors, the 
projection of range errors onto a baseline, satellite and epoch separation. A high sum 
of the satellite elevation angles for a single station can also be used to determine favor- 
able observing times, as simultaneous high elevation angles correspond to high PDOP 
Values. 

Mv results reconfirm Trimble Navigation’s [1987b] recommendation to include the 
infinite PDOP peak in the observing session; however, I could not confirm the Federal 
Geodetic Control Commuttee’s [1986] proposal to stop at a GDOP of 5.0. I could not 
identify any consistent observing stop point that improved accuracy other than to stop 
when four satellites were no longer available. 

When less than 30 minutes of data are used, the goodness of the C,A code has a 
great effect upon the carrier phase solution because C/A code positions are used to es- 
timate the carrier phase solution and to compute the time tags of the carrier phase 
measurements. The number of rejected triple difference measurements provided the best 
indicator of the quality of the C;A code estimates and the carrier phase solutions. In- 
cidences of poor accuracy of the carrier phase solution were found to be caused bv poor 
C A code estimates. The carrier phase solutions could then be improved bv correcting 
the C'A code solutions, and in turn the carrier phase time tags, by pre-processing the 
C. A code measurements or accepting more carrier phase measurements. Pre-processing 
the C/ A code measurements did not improve the accuracy of the carrier phase solutions 
when the number of rejected triple differences was less than 10% or when the solution 
standard deviation of the slope distance was less than 10 m. 

The long-term average of the best pseudorange solutions produced results conipa- 
rable to the mean carrier phase solution, except that the pseudorange solutions had 


Variances several times the magnitude of the carrier phase solutions. 


D. RECOMMENDATIONS 
The results of these case studies should be tested on survevs conducted with 


baselines of varving orientation and length and with different satellite configurations. 
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Currently, TRIMVEC does not support the use of precise ephemerides or the com- 
putation of uncorrelated triple differences. Should those capabilities be implemented in 
the future, or provided by third-party software, the data should be reprocessed and an- 
alyzed to isolate the tropospheric and ionospheric error contributions and to improve 
the baseline component results. Both of those capabilities would allow better insight 
into the eflects of geometry and the tropospheric and ionospheric errors. 

The usual role of the correlated triple difference is to aid in fixing cycle slips when 
the receiver loses lock and to provide the initial estimates for the double difference 
least-squares processing. The current data should be reprocessed in the double differ- 
ence mode to determine whether the correlated triple difference solution provides suffi- 
clent accuracy to fix the cycle slips and estimate the initial integer ambiguities in long 
baseline surveys. 

The standard values used to estimate accuracy could be improved by a more rigor- 
ous fixing of the antenna locations. Both sites should be subjected to a network ad- 
justment, either locally or simultaneously (possiblv in conjunction with the NGS VLBI 
crustal motion studies). The extra baselines required for a network solution could be 
measured using the third Naval Postgraduate School receiver. 

Studies should be conducted to determine the effects of using meteorological obser- 
vations far removed from the Sand Point antenna site. A temporary meteorological 
station could be set up near the antenna site (possibly in cooperation with the Weather 
Service Forecast Office (WSFO) located on the grounds of the Western Regional Center 
or with the nearby University of Washington). Future Naval Postgraduate School GPS 
surveys would be aided by portable meteorological instruments, preferably that made 
digital records of the temperature, humidity, and pressure. 

Because of the importance of elevation angles to accurate results, studies should be 
conducted to determine the performance of the TRIMVEC supplied tropospheric mod- 
els( modified Hopfield and Marini ) at various elevation angles. 

While this study has shown the utility of the single-station PDOP, a more complete 
DOP may provide a better satellite selection aid. Such a complete DOP ( or Dilution 
Of Relative Position (DORP)) should incorporate covariances for baseline components, 
Satellite orbital errors, receiver timing errors, ionospheric and tropospheric delays, un- 
certainties in the reference station coordinates, and cross-covariances. The DORP 


would also be formulated with respect to the type of differencing to be used. 


Lastly, the amount of data processed for these case studies was limited by the speed 
of the microcomputer and the necessity of transferring the results to the Naval Post- 
graduate School mainframe for analysis. As more GPS surveys will be conducted in the 
future, the demands for processing power will increase. It will become imperative that 
the processing software be ported to the mainframe or to a networked nuni-computer. 
Those computers will allow multiple access to the software as well as speeding the 
processing. AS an interim measure, the RAM of the current ensemble of microcomput- 
ers should be increased to several megabytes. This increase in RAM will allow the use 
of RAM-disks which will speed up the data loading - which is the most time consuming 


of the processing procedures. 
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APPENDIX. BATBLD.BAS LISTING 


10 REM PROGRAM BATBLD. BAS BY R. BOUCHARD NOV87 

20 REM PROGRAM BATBLD: BUILDS THE BATCH FILE FOR TRIMVEC PROCESSING. 

30 REM READS WX OBS FROM SEA. MET FILE 

40 REN DAY 283 REQUIRES ITS OWN BATBLD PROGRAM. 

50 INPUT "batch file name"; OFLS 

60 3۲8<" tt 

70 S1$="command /c tbf nodd. tem " 

80 OPEN OFLS FOR APPEND AS #2 

90 DIM A$(18),JD$(29),MM$(29),DD$(29),PH(29),PM(29) 

95 REM initialize Julian Day array 

100 JD$(1)="272": JDS(2)="273": JDS(3)2 274 ۰ 1:0 (2 

: JDS(5)="276" 

110 JD$(6)="279": JDS(7)="280": JD$(8)="281": JDS(9)="282" 

: JDS(10)="287" 

120 JD$(11)="290": JD$(12)="293": JD$(€13)="294": JD$(14)="295" 
JD$(15)2'297" 

130 JD$(16)2' 300"; JDSCIZ)Z 301 | 0ص 9 0 خر ۱9 اولال‎ ٦ 

140 JD$(20)="308": JD$(21)='309": JDS$(22)='314" 

145 JDS(23)="315": JDS(24)="317": JDS(25)="318” 

147 JD$(26)="325": JD$(27)="329" : JDS(28)="321" 

150 MMS(1)="09": MMS(2)="09" 

160 FOR I=3 TO 18: MMS(I)="10": NEXT I 

170 FOR 1219 TO 28: MMS$(I)="11": NEXT I 

180 DD$(1)="29": DDs(2)="30": DDSs(3)="01": DDs(4)="02": DDs(5)="03" 

: DD$(6)2'06": DDS( 

7)="07": DDs$(8)=""08": DD§(9)="09": DD$(10)="14": DD$(11)="17" 
DDs(12)="20" 


190 DD$(13)= 21": DD$(14)= 22: DD$(15)="24': DD$(16)= "27": DD$(17)="28" 
200 DD$(18)="29"': DDS(19)="03": DD$(20)="04": DDS(21)="05": DDS(22)="10" 


205 DD$(23)="11": DD$(24)="13": DD$(25)="14" 
207 DD$(26)="21": DDS(27)2" 25": DDs(28)="17" 
210 AS(1)='a01": A$(2)="a02": AS$(3)="a03": AS(4)="a04": AS(5)="a05" 


: A$(6)2'a06" 
220 AS(7)="a07'": A$(8)="a08": AS(9)="a09": AS(10)="a10": AS(11)="a11" 
As(12)="a12" 


230 As(13)="a13": AsS(14)="a14": AS(15)="a15": A$(16)="a16" 
AS(17)="a17'": As(18)="a18" 

240 FOR I=1 TO 7: PH(I)=19: NEXT I 

250 FOR I=8 TO 13: PH(1)=18: NEXT I 

260 FOR I=14 TO 21: PH(1)=17: NEXT I 

270 FOR I=22 TO 25: PH(I)=16: NEXT I 

275 FOR I=26 TO 28: PH(I)=15: NEXT I 

280 PM(1)=34: PM(2)=31: PM(3)=26: PM(4)=23: PM(5)=18: PM(6)=5: PM(7)=1 

290 PM(8)=58: PM(9)=54: PM(10)=33: PM(11)=22: PM(12)=8: PM(13)=3 
PM(14)=59 

300 PM(15)=51: PM(16)=39: PM(17)=35: PM(18)=31: PM(19)=14 

310 PM(20)=6: PM(21)=2: PM(22)=42 

315 PM(23)=37: PM(24)=29: PM(25)=24: PM(26)=56: PM(27)=39: PM(28)=8 

317 REM begin building the batch file for each day. 

320 FOR IK=1 TO 28 

330 OPEN "di 1, sea. met: 
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INFS="sa"+JDS(IK)+" ma"+JDS( IK) 

S2S=S1S+INFS 

PRINT EOF JA... y JDSCIK) 

INPUT "start hour"; SH 

INPUT “start minute"; SM 

۶۰۷۰۰٠۳۰ 105۶۱۷۷ LK) PMCIK) 

SHS=STRS$ (SH) 

SMS=STR$(SM) 

۲6-0 

1-1 

REM Find the WX OB 

Mei ees these iC 1). T1C1)sRICL) ,P2CT),T201),R2(1) 
IF(CIDSZJDSCIK))AND(SH-VAL(IHS))) THEN 470 

GOTO 440 

ET ZEIT 2D E T2(I)s(TI2(I)-32)5579 

I=I+1 

( 1 )۴۱2 ر1 )12( 1 )۶2 ار 10و( 1 )لل و(1 :11۹5,115 ۸1 
THEN 0‏ ( ر 1۳ )لالج ۱۱۳۱۲ 

KEN Subtract 2 mb from SEA ob 

P1(1)=P1(1)-2! 

REM Convert NPS Temp to Celsius 

0ٰ 029 

GOTO 480 

IE=I-1 

REM interpolate WX OBS TO GPS START TIME 

PETOTOS PICIVECPIC2YSPRI(I))*PC 
EXIME LM ٥ PC 

لل( (1 )1 2 ۲۸۱۸ ۱۱۱ ۷۲۱ 1 

MASIA ۱1 2 ۶ ۱) 1 (۵ 
RR1(1)=R1(1)+(R1(2)-R1(1))*PC 

AI رز 1 )2 (2 )2 ل1‎ PC 

REM Compute running mean of first hour 

BEA E ) ۳۳۰۸ ۳ 122 

EE PRAG 72 

۱۳۱۱ 2 (5۱1) ۳۳۱2۵2 

۷۴۰۰٠٠٠12 د‎ 

RR1(2)=(RR1(1)+R1(2))/2 

2 2 گ۰ 2 ۱2۱ 

WT=1-PC 

REM Compute remaining running means 

POR IS 3 TO IE 
IA=I-1 

IWT=1+WT 
PRICI)=(PRICIA)*WT 
۳۲2 ) 1 (<) ۳۶2 ) 1۸ (۷ 
1۳1) 1 (<6 1۳1) 1۸ (1 
1۳۴2) 1 (2) 1۳52 ) 1۸ ( 
11۳1 ) 1 (2) ۳1) 1۸ (1 
RR2(I)SC(RR2(IA)*WT 
WI-IWT 

NEXT I 

M1=SM 

H=SH 

C=0 

I-2 

DHRZPH(IK)-SH: DMIN-PM(IK)-SM 


) ۱) 1۵( ٤٣۱) 1(((/2( INT) 
((P2(IA)+P2(1))/2))/CIWT) 
(TALA LO AT) 
AAA 7-2 oy GWT) 
(CRICIA)+R1IC1I))/2))/CIWT) 
((R2CIA)+R2(I))/2))/(CIWT) 


4+ + + + + + 
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340 
350 
360 
31710 
380 
29 
400 
410 
420 
430 
435 
44 

450 
460 
470 
480 
490 
500 
505 
510 
DS 
520 
530 
540 
550 
560 
570 
580 


IF(DMIN<0) THEN DHR=DHR-1: DMIN=DMIN+60 


NP=DHR*6+(DMIN/10) 


REM Interpolate hourly running means to 10 minute intervals 


I=I+1 


I=I+1 


FOR IL=1 TO NP 
M1=M1+10 

IF(M1>60) THEN M1=M1-60: H=H+1: 
WT=M1/60 

IF(M1=60) THEN M1=M1-60: H=H+1: 
MP=M1 
HP=H 

IWT=1-WT 

TA=I-1 

۲۲۷ 12۳161 ) ۲۸ ( ۰۳) 1] (+۳1) 1۷ 
PW1=INT((10*PW1)+.5)/10 

۳۷۷ 22۳1۳2) 1۸ (۰۳) 1۷] ( +۳۴۲2) 1 
PW2=INT((10*PW2)+.5)/10 
TWIZTRICIA)*(IWT)TTRI(I)*WT 
TW1=INT((10*TW1)+.5)/10 
TW2=TR2( IA)*( IWT)+TR2(I)*WT 
TW2=INT((10*TW2)+.5)/10 

1۷ 12۳۳1) 1۸ (۰۰6 1۷1 ( +۳۳1) (۷۲ 
RW1=INT((10*RW1)+.5)/10 
RW22RR2(IA)*(IWT)TRR2(I)*WT 
RW2=INT((10*RW2)+.5)/10 

REM convert to strings 
PS1$5=LEFTS(STRS(PW1),7) 
PS2S=LEFTS(STRS(PW2),7) 
15152101195 315) ال‎ 5) 
TS2S=LEFTS(STRS(TW2),5) 
RS15=LEFTS(STRS(RW1),5) 

۲9 2 9 21 ۱ 518 5) ٨٩١۷ 2( 5) 
HSS=STRS(HP) 
S35=S25+SPS$+AS(IL)+". "+JDS( IK) 
MSS=STRS (MP) 


PRINT $3$+PS1$+TS1$+RS1$+PS2$+TS25+RS2$+SPS+MM$(IK)+SP$ 


) 1٨۹۳ ۶ 


850 
860 
865 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
290 
1000 
1010 
1920 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1159 
1140 
115% 
1160 
1170 
1180 
+DD$ 


1190 PRINT#2, S3$+PS1$+TS15+RS1$+PS2$+TS2$+RS2$+SP$+MM$(IK) 


+DDS( IK)+HSS+MSS+SHS+SMS 


TOES 


1200 LPRINT S$3$+PS1$+TS1$+RS1$+PS2$+TS2$+RS2$+52$+MM$(IK)+SP$ 


56 


(IK)+HSS+MSS+SHS+SMS 
NEXT IL 
CLOSE #1 
NEXT IK 


+005 
1210 
1220 
1230 


1240 CLOSE #2 


END 
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